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Abstract 


Creation of Cold Dark Matter (CCDM), in the context of Einstein Field Equations, leads 
to a negative creation pressure, which can be used to explain the accelerated expansion of the 
Universe. Recently, it has been shown that the dynamics of expansion of such models can not 
be distinguished from the concordance ACDM model, even at higher orders in the evolution of 
density perturbations, leading at the so called “dark degeneracy”. However, depending on the 
form of the CDM creation rate, the inclusion of spatial curvature leads to a different behavior of 
CCDM when compared to ACDM, even at background level. With a simple form for the creation 
rate, namely, T oc we show that this model can be distinguished from ACDM, provided the 
Universe has some amount of spatial curvature. Observationally, however, the current limits on 
spatial flatness from CMB indicate that neither of the models are significantly favored against 
the other by current data, at least in the background level. 
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I. INTRODUCTION 


Since 1998, the series of observations of the luminositv-redshift relation of the Super¬ 
novae Type la (SNe la) has set a new era in cosmology |l|-|4j]. Those observations, com¬ 
plemented by the observations of the Cosmic Microwave Background (CMB) anisotropies, 
Baryon Acoustic Oscillations (BAO), Hubble parameter in different redshifts, strongly sug¬ 
gest that the Universe has a great order of spatial flatness and has entered in a late phase 
of accelerating expansion |3|-|8|. 

Inside the context of relativistic cosmology, the accelerated expansion of the Universe 
is usually attributed to a new dark component called dark energy which possesses as 
main feature the negative pressure. The most favoured candidate for dark energy is the 
cosmological constant, A. This new component not only can fit the SN la observations 
but can also recover the flatness of the spatial hypersection of the Universe, as predicted 
by inflationary theory and corroborated by CMB observations. 

Nevertheless, this concordance model has some serious issues. For example, the density 
of cosmological constant A, responsible for the late acceleration of expansion of the Universe 
as estimated from observational data (p\) is almost null (10“ 47 g/ern 3 ) and must be fine 
tuned to explain quantitatively the acceleration. Due to its equation of state = — Pa), 
the cosmological constant could rise from the vacuum energy of quantum fields. However, 
the theoretical estimation to the vacuum energy through quantum held theory ( p v ,Th ) is 
up to 122 orders of magnitude bigger than the observational density [9]. 

Inside the General Relativity, the negative pressure is the key element for acceleration. 
This can occur naturally in thermodynamical process departing from the equilibrium, for 
instance, the matter creation process at expenses of gravity. This phenomenon gives rise to 
a term of negative pressure which should be considered at the level of Einstein Equations, 
as shown by Prigogine and collaborators and formulated in a manifestly covariant way by 
Lima, Calvao & Waga |10l. 111! ]. The inclusion of the backreaction at the level of Einstein 
Field Equations was determinant to the rise of a new class of cosmological models with 
matter creation. 

Many CCDM models have been proposed in the literature, each of those have different 
phenomenological creation rate dependencies. Among them, recently, a model from Lima, 
Jesus and Oliveira (LJO) was interesting because it was shown that this model leads to the 
same background evolution as the concordance ACDM model, for any spatial curvature. 
Further, it was shown that even at linear density perturbation evolution, such a degeneracy 
persisted, leading to the so called “dark degeneracy”. 

A first alternative to break such a degeneracy was given by Jesus and Pereira 12], which 
have found, directly from quantum held calculations, a creation rate which depended on 
the dark matter mass. They have used observational data to constrain the dark matter 
mass and have argued that with better constraints in the future, this could be used to 
distinguish CCDM from ACDM. 
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In this paper we investigate a cosmological model driven only by cold dark matter 
(CDM) creation at expenses of gravitational field in which the rate of CDM creation 
evolves reciprocally to the expansion rate, and we include the possibility of nonzero spatial 
curvature. We assume that the created particles are described by a real scalar field and 
consequently the created particles are its own antiparticles. Similarly to the standard 
model, the scenario presented here has the same degree of freedom and is also capable 
of explaining the accelerating expansion. We show that, with some amount of spatial 
curvature, the CCDM behavior differs from ACDM and both can be distinguished using 
observational data. 

In Section II, we discuss the dynamics of the universe with the pressure due to creation. 
In Section III, we discuss an specific rate of dark matter creation. In Section IV, we 
constrain the free parameters of the model. In Section V, we compare the model with 
other models on the literature. Finally, we summarize the main results in conclusion. 


II. COSMIC DYNAMICS ON MODELS WITH CREATION OF CDM PARTI¬ 
CLES 

We will start by considering the homogeneous and isotropic FRW line element (with 
c = 1): 

ds 2 = dt 2 — a 2 (t ) (— 

\ 1 — kr z 

where k can assume values —1, +1 or 0. 

In this background, the Einstein field equations are given by 

a 2 k 

8irG(p ra d + Pb + Pdm ) = 3—^ + 3—, (2) 

er a z 

and 

8nG(p rad + p c ) = -2- K- (3) 

a a- a z 

where p rad , Pb and p d m are the density parameters of radiation, baryons and dark matter, 
Prad = Prad /3 is the radiation pressure and p c is the creation pressure. 

For the radiation and baryon components, the energy conservation are given by the 
usual expressions: 

a 

Prad "F 4 Prad 0, (4j 

a 

and 

p b + 3-pb = 0. (5) 

a 


+ r 2 dd 2 + r 2 s'm 2 6dq F 


( 1 ) 
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where each overdot means one time derivative and we have used that p rac t = p ra d /3 and 
Pb = o. 

On the other hand, when the creation process is considered_we should take into account 
a matter creation source at level of Einstein Field Equations 
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( 6 ) 


where T is the rate of dark matter creation in units of (time) -1 . 

As shown by H[, the creation of dark matter is responsible for an extra pressure term 
at the level of Einstein Equations, the so called creation pressure, p c , 


Pc = 


PdmX' 

’ 3 H ’ 


( 7 ) 


where we have considered an “adiabatic” creation, i.e., the case when the entropy per 
particle is constant. 

As a consequence of the above equation, one can see that the dynamics of the universe 
is directly affected by the rate of creation of cold dark matter, T. In particular, in the case 
T > 0 (creation of particles) we have a negative pressure creation and in the case r —> 0 
we recover the well known dynamics when the universe is lately dominated by pressureless 
matter (baryons plus dark matter). 


III. CREATION OF COLD DARK MATTER (CCDM) MODEL 

The difficulties in identifying the nature of dark energy led the cosmologists to a quest 
for better candidates to explain the late acceleration of the Universe. In the literature, 
models with CDM creation has been discussed as a viable explanation to this recent phe¬ 
nomenon. It has been shown that under a convenient choice of the particle creation rate 
T, this scenario is able to support the observed dynamics, linear structure formation and 
thermodynamics features of the late Universe jTl . 
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We argue that a natural dependence of the CDM creation rate is on the expansion rate, 
as already proposed in other CCDM models present on the literature. It has already been 
studied a linear dependence with the Hubble parameter [l5j, and at the time of writing 
this paper, it has been proposed a power law dependence on Hubble parameter |16|. It is 
clear to us that the expansion acceleration is a recent feature, so the CDM creation also 
must be recent. As the Universe evolves, the Hubble parameter decreases, and T must 
increase. We can satisfy those conditions by assuming T to be a negative power law of H, 
or, in the simplest case, T oc if -1 . So, in this work, we consider the following creation 
rate: 


r 


m 

3 “ir 


( 8 ) 
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where a is a constant free parameter of the model which drives the creation rate and the 
factor 3H 2 has been introduced for mathematical convenience. This is also interesting 
because, as we shall see, with the identification a = 12a, the flat ACDM is a particular 
case of our CCDM model, when we neglect the baryon contribution. 

Since we are considering only the late phase of the dynamics of the universe, we can 
neglect the radiation terms from now on. Thus, by combining Eqs. (J2J) and (|3]) , we have 


a 

a 


4vr G 


(Pb T Pdm T 3 Pc) 


Replacing p c from Eq. (J7J) , we may write 


a 

a 


4ttG 

^T 


Pb T Pdm 



(9) 


( 10 ) 


Using that - = H + H 2 and changing variables from time to redshift, we End 

dH _ H H 2 Q b (l + z) 2 T H 2 -H 2 Q k {l + z) 2 / T\ 
dz 1 + z 2 H 2 2H(1 + 0 ) V Hj 


where we have used the solution of ([5]) to baryon density, p b = p b o(l + z) 3 , h2 b = is 
the present baryon density parameter, and Vt k = is the present curvature density 
parameter. Replacing the creation rate (jHJ) and changing to dimensionless variable y = 

(f) 2 , we find 


dy_ 

dz 


0 y - a t 3afh,(l + Z Y 

o- H 
1 + z y 


+ z) 



( 12 ) 


which can not be solved analytically, except for some special cases. For example, if the 
Universe is spatially flat and we can neglect baryons = 0), we can solve (112)1 to 

End: 

H{z) 2 = H*[a + {l-a)(l + z) 3 } (13) 


which corresponds to the flat CCDM LJO model, or to flat ACDM model, with the identi¬ 
fication a = Vt a. However, by introducing baryons or curvature, this model can not recover 
the LJO or the ACDM model anymore. So, this model can be useful to discriminate be¬ 
tween CCDM models and ACDM model even in the background level. On the other hand, 


the LJO and ACDM models can only be distinguished at the pertur 


oation levels 


only in the absence of the separation of dark matter components 17]. 


and 


There is also one more analytical solution. If we neglect curvature (14 = 0), retaining 
baryon density (14 ^ 0), we find an analytical solution to (1T2]) : 


rr \ 2 

w j =m + z ) 3 


1 + w 


(1 - n b )e n » 
14 


-i' 


(14) 
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FIG. 1: Numerical solutions for H/Hq in function of the redshift 0 , and its sensitivity to the free 
parameters a and Left) Evolution of H(z) with a = 0.7 and = 0 (solid curve), = 0.6 
(dashed curve) and fl*, = —0.5 (dashed-dotted curve). Right) Evolution of H(z) with fl*, = 0 
and a = 0.0 (solid curve), a = 0.7 (dashed-dotted curve) and a = 0.9 (dashed curve). 

where W ( x ) is the principal Lambert W function, also known as product logarithm, real 
solution to equation x = W(x)e w ^ x \ 

However, if we neglect baryons only (Of, = 0, Q*, ^ 0), or if we consider the full equation 
no, with Of, ^ 0 and O*. ^ 0, we can not find an analytical solution to H(z), and we have 
to resort to numerical methods. Nevertheless, if curvature and baryonic contributions can 
both be considered small (0 < 0& 1, |Ofc| <C 1), we can find an approximation, 

1 = « + (1 — a )(l + z ) 3 + — (1 + z ) 2 ] + 30*, log(l + z) + 

+ f (j~~^ ~ lQ g [ a + ( 1 ~ a )( 1 + ^) 3 ] (15) 

which we have found by solving (TT2j) with Q b — Q k — 0, replacing the solution on (1121) and 
solving again, with Of, ^ 0, O fc ^ 0. 

In Figure [D we show the numerical solutions of H(z)/Hq for some values of the free 
parameters of the model, namely, 0& and a. It is worthy to remark that the model has as 
particular cases two well known models: the Einstein-de Sitter, for O fc = 0 and a = 0 and 
the flat ACDM for O fc = 0 and a ~ 0.7 and Of, = 0. 


IV. OBSERVATIONAL CONSTRAINTS 


In this section, we obtain constraints to the free parameters of the model, namely, 0& 
and a. In order to do this, we considered some measurements of the Hubble parameter, 
H(z) [18] and the Supernovae Type la (SN la) dataset of Union 2.1 
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A. H{z) Constraints 


Hubble parameter data as function of redshift yields one of the most straightforward 
constraints on cosmological models. In general, these observational data depend on astro- 
physical assumptions, not depending on any background cosmological model. However, is 
not easy to obtain such data, as it is not directly observable, but is rather inferred from 
astrophysical observations. 

At the present time, the most important methods for obtaining H(z) data are (i) through 
“cosmic chronometers”, for example, the differential age of galaxies (DAG), (ii) measure¬ 
ments of peaks of acoustic oscillations of baryons (BAO) and (iii) through correlation 
function of luminous red galaxies (LRG) [8|. In this work, we use the data compilation 
of H(z) from Farooq and Ratra [l8|], which is, currently, the most complete compilation, 
with 28 measurements. 

From these data, we perform a y 2 -statistics, generating the Xh function of free param¬ 
eters: 

H 0 E(zi, cr, Dfc, f£) Hi 


Xh = 


28 

£ 

i=l L 


u Hi 


(16) 


where E(z) = and H(z) is obtained by solving numerically Eq. (I12[) . 

Throughout every analysis on this paper, we fix the baryon density parameter at the 
value estimated by Planck and WMAP: = 0.049 |6 h a value which is in agreement with 

Big Bang Nucleosynthesis (BBN), as shown on Ref. [191 ]. 

As the function to be fitted, H(z) = HqE(z), is linear on the Hubble constant, Hq, we 
may analitically project over H 0 , yielding Xh : 


Xh = C- 


B 2 

~A 


where A = £" = i ^ B = ELi C = ELi yj and E i = 
analysis can be seen on Figure [2] (left). As can be seen, the results from H(z) data alone 
yield very loose constraints on the plane VLk-ot. In fact, over the region 0 < a < 1.4 
and —1 < flfc < 1, only the 68% c.l. statistical contour could close. The minimum y 2 was 
Xmin = 16.269, which is too low for 28 data, yielding a y 2 per degree of freedom y 2 = 0.626. 
The best fit parameters were a = 0.79llo;gg 5 , = in the joint analysis. We 

believe that such a loose constraint can be due to a underestimate of uncertainties on the 
H(z ) compilation data, which is evidenced by its low y 2 . This issue has been addressed by 


H? 


— H{zj) 

H 0 ' 


(17) 


The result of such 


181 ] by combining H{z) data with other constraints. While 18] uses a prior over H 0 , we 


choose to use a prior over Qk, as a prior over Hq did not affect much the constraints found 


with H 
results 


z) data only. We have considered a prior over flk based on Planck and WMAP 

A 


Planck + WMAP indicate = —0.037_ ao49 , at 95% c.l., in the context o 


ACDM. 

Based on this result, along with the symmetrization process suggested by 20j], we use 
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a prior of Q*. ± a U}k = —0.043 ± 0.046. We will refer to this simply as the CMB prior. 
The results can be seen on Figure [2] (right). As can be seen there, the limits over a 
and Q k are quite better. We have found X 2 mxn = 16.344, a = 0.775lg;o64laii 8 iai6 an d 
Q k = —0.041 ± 0.069 ± 0.11 ± 0.16. In order to improve these constraints, we made a 
combined analysis with SN la data. 


B. Supernovae Type la Bounds 


The parameters dependent distance modulus for a supernova at the redshift z can be 
computed through the expression 


/x(z|s) = m — M = 5 \ogd L + 25, 


(18) 


where m and M are respectively the apparent and absolute magnitudes, s = ( H 0 , ct, Q k ) 
is the set of the free parameters of the model and d k is the luminosity distance in unit of 
Megaparsecs. 

Since in the general case H(z) has not an analytic expression, we must define d,L through 
a differential equation. The luminosity distance g(l can be written in terms of a dimen¬ 
sionless comoving distance D by: 


d L = (l + z)^D 
c 


(19) 


The comoving distance can be related to H(z), taking into account spatial curvature, 
by the following relation 21[: 


(H\ 2 _ £l k D 2 + 1 
\H 0 ) ~ V ~ D' 2 


( 20 ) 


where the prime denotes derivation with respect to redshift £. Inserting this relation in 
Eq. (fT2lh we obtain a differential equation for D\ 

= n k DD' 2 _ 3 D’ D /3 

1 + fl k D 2 2(1 + z) + 2(1 + fl k D 2 ) 

3afi & (l + z) 2 D' 5 
2(1 + Vl k D 2 ) 2 ’ 

To solve numerically this equation we have used as initial conditions 71(0) = 0 and 
-D'(O) = 1. In order to constrain the free parameters of the model we considered the Union 
2.1 SN la dataset from Suzuki et al. [s|. The best-fit set of parameters s was estimated 
from a y 2 statistics with 


3ck 
1 + z 


+ flfc(l + z) ( 1 


3 aD 
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1 + fl k D 2 















FIG. 2: The results of our statistical analysis, with contours for 68.3%, 95.4% and 99.7% confi¬ 
dence intervals. In both panels, the dashed lines correspond to SNs Union 2.1 data. Left) Dotted 
lines: constraints from H(z) data, solid lines: constraints from SNs + H(z). Right) Solid lines: 
constraints from SNs + H(z) + CMB prior on 0*,. 


where n % (z |s) is given by (IT8j) . /%(z) is the corrected distance modulus for a given SNe la 
at Zi being cq its corresponding individual uncertainty and N = 580 for the Union 2.1 data 
compilation. 

In Figure [2] (left), we display the space of parameters Qk~ ol and the contours for la, 2a 
and 3cr of confidence intervals. Using SN data alone, and marginalizing the nuisance pa¬ 
rameter h (Ho = lOOh km s^ 1 Mpc _1 ), we constrain the free parameters as a = 0.7881qq§ 6 , 
hi*, = 0M±° o il at 68% confidence level and xl — 0.974. We also show, on Figure [2] (left), 
the combination SNs la + H(z), obtained by adding Xs N + Xh an d marginalizing over 
h. As can be seen, adding H(z) data alone to SN la data, barely changes the SN con¬ 
straints. In fact, we have found, on this case, a = O.792+q;q0 1 +q;q9 3 , Xl k = ^-^-o'II-oao, 
with xl — 0.956. The main difference on this case, is that the 95% c.l. contours could 
close on the region considered. 

So, in order to find the best possible constraints with the data available, we made the 
full combination of SNs + H(z) + CMB prior on fi*,. In this case, the constraints were 
quite restrictive, as shown in FigureE] (right). As one can see, the CMB prior makes a great 
cut over the Union 2.1 contours, as is strongly constrained, in this case. We have found, 
on this case, in the joint analysis: a = 0.768lo:o3mao5mao69, ^k = -0.036lo:o68-aii-ai5> 
xl = 0.957. Table H] summarizes the analysis results, including the analysis with the Ho 
prior, based on the current limit on Ho given by 22], Hq = 72.0 ± 3.0 km s~ x Mpc 


-l 
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Data 

a 

^k 

xl 

H(z) 

0 791+ 0 - 18 
u ' (yi -0.085 

n 04+ 0 - 46 
u - u ^-0.40 

0.626 

SNs 

n 7 qq+0-23 

U./55 o 086 

0 04+ 0 - 39 
u - u ^-0.35 

0.973 

SNs +H(z) 

Pi 7QO+0-H +0.26 
u ‘ * y —0.061—0.093 

n n f v _ i _ 0-25+o.42 
U - UO -0.24-0.40 

0.955 

SNs +H(z) + Hq prior 

O 707 +0.11 +0.25 
u *' 0/ -0.065-0.096 

n rn +0.26+0.43 
u * ul -0.26-0.42 

0.957 

H(z ) ± Hq prior 

f) 7Qft+0-17+0.61 

U. /OU o 085-0.12 

n 

u * uu -0.35-0.57 

0.649 

H(z)+ CMB prior 

rv 'y'yr+0.059+0.098+0.14 
u - mo -0.064-0.11 -0.17 

-0.041 ±0.069 ±0.11 ±0.16 

0.629 

SNs+i7(,2)+CMB prior 

n 7£io+0.034+0.058+0.084 
U. /DO_o 031-0.051-0.069 

n n+o.15 

U.UoO_ 0 068—0 .11— 0.15 

0.955 


TABLE I: Results of the joint analysis for the different combinations of data. Limits on the 
parameters correspond to 68.3%, 95.4% and 99.7% c.l. as explained on text. 


V. COMPARISON WITH CCDM AND ACDM MODELS 


In the absence of baryons and if the spatial curvature vanishes, this model coincides both 
with the concordance ACDM model and LJO model 131 ]. where r oc + Also, this model 
has the same creation rate as the CCDM1 model of Ref. 23]. However, in their analysis, 
they have not included baryons nor considered nonzero spatial curvature. As baryons 
are separately conserved, the DM creation rate T oc A leads to a huge difference on the 
Universe evolution, as we have shown on our analytical solution, Eq. (fT4lb Numerically, 
taking baryons into account, for ~ 0.04, leads to a relative difference of about 10% on 
H(z), for z>l. 

It has already been shown that the LJO model can not be distinguished from ACDM, 


for any value of spatial curvature, neither at background level 13] nor at perturbation 


r-L 


level [17j. In this manner, LJO gives rise to the so called “dark degeneracy” 17], where, 


through cosmological observations, one can not determine if it is the quantum vacuum 
energy contribution (ACDM) or the quantum vacuum dark matter creation (CCDM) which 
accelerates the Universe. However, in the model proposed here, if the Universe has some 
amount of spatial curvature, we can distinguish it from ACDM. 

In fact, we have proposed a natural dependence of the creation rate over the expansion 
rate, and the direct inclusion of spatial curvature already breaks the “dark degeneracy”. 
The question now is: can the LIniverse be nonflat enough in order to distinguish both 
models? Part of the answer is on Figure [2] As we have seen, SNs +H(z) alone, which 
are observational data quite independent from cosmological models, are not enough for 
constraining the spatial curvature. Namely, looking at its 95% confidence limits over the 
curvature, we have —0.19 < < 0.30 so, the Universe can be flat, quite open (D^ ~ 0.30) 

or quite closed (LR. ~ —0.19). Spatial curvature on the border of this limit is enough for 
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FIG. 3: Comparison between H(z) curves from CCDM and ACDM for some values of 0^. The 
regions correspond to 68% c.l. deviation over the Hq value, as explained on the text. Also shown 
is the H (z) data used for the statistical analysis. 


distinguishing CCDM from ACDM, as can be seen on Figure [3] 

We have plotted, on Figure El H(z) curves for some values of D*,, in the context of 
CCDM and ACDM. In all curves, we have fixed a = Da = 0.792, the CCDM best fit from 
SNs +H(z), and we have used the H 0 value from 22], H 0 = 72.0 ± 3.0 kms~ 1 Mpc" 1 , with 
the regions corresponding to 1 -a variation over the best fit. As we can see, if the Universe 
is closed, CCDM and ACDM are harder to be distinguished, with a small distinction 
only at high redshifts. If, however, the Universe is open, the distinction can be made at 
intermediate redshifts, and, at high redshifts, the distinction is clear even with the current 
set of H(z) data. 

However, if we take into account the prior over given by CMB, we can not distin¬ 
guish them, as we have, in this case, an strict 95% confidence limit over the curvature of 
—0.104 < < 0.030. This certainly is not enough to distinguish CCDM from ACDM. 

Here, one could think that CMB constraints on ACDM could not be used, directly, to 
constrain CCDM models, as the former model has not matter creation and the latter has. 
However, for the creation rate used here, T oc -A, creation is negligible on early Universe 
evolution, thus not changing the signatures imprinted on the last scattering surface. Thus, 
we conclude that the Universe is quite spatially flat, even in the context of this partic¬ 
ular CCDM model, so we can not distinguish it from ACDM. The “dark degeneracy” is 
maintained from this analysis. 
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VI. CONCLUSION 


In this work we proposed a new cosmological model, based on the matter creation 
phenomena. We have shown that the proposed model is able to explain the background 
accelerating dynamics of the Universe, without a new dark fluid with negative pressure. 

We have demonstrated that the present model is able to avoid the “dark degeneracy” 
simply through the presence of a baryonic content or the spatial curvature. For both cases, 
the model presents a distinguishable Hubble expansion from the ACDM. 

We have shown that this model can be distinguished from LJO with the inclusion of 
baryons and spatial curvature, including a relative difference of about 10% on H(z), for 
z > 1 and U;, ~ 0.04, as mentioned on Section V. Nevertheless, the current data from SNe 
la and H(z) are still not able to observationally favor any of both models in the background 
level, as we have discussed. 

Further investigations including perturbation analysis should be made, in order to re¬ 
alize if more significant distinctions between CCDM and ACDM models can be found at 
higher orders of density perturbations evolution. 
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